Identifying Candidate Gene Drivers Associated with Relapse in Pediatric T-Cell Acute Lymphoblastic Leukemia Using a Gene Co-Expression Network Approach.

Pediatric T-cell Acute Lymphoblastic Leukemia (T-ALL) relapses are still associated with a dismal outcome, justifying the search for new therapeutic targets and relapse biomarkers. Using single-cell RNA sequencing (scRNAseq) data from three paired samples of pediatric T-ALL at diagnosis and relapse, we first conducted a high-dimensional weighted gene co-expression network analysis (hdWGCNA). This analysis highlighted several gene co-expression networks (GCNs) and identified relapse-associated hub genes, which are considered potential driver genes. Shared relapse-expressed genes were found to be related to antigen presentation (HLA, B2M), cytoskeleton remodeling (TUBB, TUBA1B), translation (ribosomal proteins, EIF1, EEF1B2), immune responses (MIF, EMP3), stress responses (UBC, HSP90AB1/AA1), metabolism (FTH1, NME1/2, ARCL4C), and transcriptional remodeling (NF-κB family genes, FOS-JUN, KLF2, or KLF6). We then utilized sparse partial least squares discriminant analysis to select from a pool of 481 unique leukemic hub genes, which are the genes most discriminant between diagnosis and relapse states (comprising 44, 35, and 31 genes, respectively, for each patient). Applying a Cox regression method to these patient-specific genes, along with transcriptomic and clinical data from the TARGET-ALL AALL0434 cohort, we generated three model gene signatures that efficiently identified relapsed patients within the cohort. Overall, our approach identified new potential relapse-associated genes and proposed three model gene signatures associated with lower survival rates for high-score patients.

Recessively Inherited Deficiency of Secreted WFDC2 (HE4) Causes Nasal Polyposis and Bronchiectasis.

RATIONALE: Bronchiectasis is a pathological dilatation of the bronchi in the respiratory airways associated with environmental or genetic causes (e.g., cystic fibrosis, primary ciliary dyskinesia and primary immunodeficiency disorders), but most cases remain idiopathic. OBJECTIVES: To identify novel genetic defects in unsolved cases of bronchiectasis presenting with severe rhinosinusitis, nasal polyposis, and pulmonary Pseudomonas aeruginosa infection. METHODS: DNA was analyzed by next-generation or targeted Sanger sequencing. RNA was analyzed by quantitative PCR and single-cell RNA sequencing. Patient-derived, cells, cell cultures and secretions (mucus, saliva, seminal fluid) were analyzed by Western blotting and immunofluorescence microscopy, and mucociliary activity was measured. Blood serum was analyzed by electrochemiluminescence immunoassay. Protein structure and proteomic analyses were used to assess the impact of a disease-causing founder variant. MEASUREMENTS AND MAIN RESULTS: We identified bi-allelic pathogenic variants in WFDC2 in 11 individuals from 10 unrelated families originating from the United States, Europe, Asia, and Africa. Expression of WFDC2 was detected predominantly in secretory cells of control airway epithelium and also in submucosal glands. We demonstrate that WFDC2 is below the limit of detection in blood serum and hardly detectable in samples of saliva, seminal fluid, and airway surface liquid from WFDC2-deficient individuals. Computer simulations and deglycosylation assays indicate that the disease-causing founder variant p.Cys49Arg structurally hampers glycosylation and thus secretion of mature WFDC2. CONCLUSIONS: WFDC2 dysfunction defines a novel molecular etiology of bronchiectasis characterized by the deficiency of a secreted component of the airways. A commercially available blood test combined with genetic testing allows its diagnosis. This article is open access and distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/).

Detecting subtle transcriptomic perturbations induced by lncRNAs knock-down in single-cell CRISPRi screening using a new sparse supervised autoencoder neural network.

Single-cell CRISPR-based transcriptome screens are potent genetic tools for concomitantly assessing the expression profiles of cells targeted by a set of guides RNA (gRNA), and inferring target gene functions from the observed perturbations. However, due to various limitations, this approach lacks sensitivity in detecting weak perturbations and is essentially reliable when studying master regulators such as transcription factors. To overcome the challenge of detecting subtle gRNA induced transcriptomic perturbations and classifying the most responsive cells, we developed a new supervised autoencoder neural network method. Our Sparse supervised autoencoder (SSAE) neural network provides selection of both relevant features (genes) and actual perturbed cells. We applied this method on an in-house single-cell CRISPR-interference-based (CRISPRi) transcriptome screening (CROP-Seq) focusing on a subset of long non-coding RNAs (lncRNAs) regulated by hypoxia, a condition that promote tumor aggressiveness and drug resistance, in the context of lung adenocarcinoma (LUAD). The CROP-seq library of validated gRNA against a subset of lncRNAs and, as positive controls, HIF1A and HIF2A, the 2 main transcription factors of the hypoxic response, was transduced in A549 LUAD cells cultured in normoxia or exposed to hypoxic conditions during 3, 6 or 24 h. We first validated the SSAE approach on HIF1A and HIF2 by confirming the specific effect of their knock-down during the temporal switch of the hypoxic response. Next, the SSAE method was able to detect stable short hypoxia-dependent transcriptomic signatures induced by the knock-down of some lncRNAs candidates, outperforming previously published machine learning approaches. This proof of concept demonstrates the relevance of the SSAE approach for deciphering weak perturbations in single-cell transcriptomic data readout as part of CRISPR-based screening.

BABA-induced pathogen resistance: a multi-omics analysis of the tomato response reveals a hyper-receptive status involving ethylene.

Prior exposure to microbial-associated molecular patterns or specific chemical compounds can promote plants into a primed state with stronger defence responses. ß-aminobutyric acid (BABA) is an endogenous stress metabolite that induces resistance protecting various plants towards diverse stresses. In this study, by integrating BABA-induced changes in selected metabolites with transcriptome and proteome data, we generated a global map of the molecular processes operating in BABA-induced resistance (BABA-IR) in tomato. BABA significantly restricts the growth of the pathogens Oidium neolycopersici and Phytophthora parasitica but not Botrytis cinerea. A cluster analysis of the upregulated processes showed that BABA acts mainly as a stress factor in tomato. The main factor distinguishing BABA-IR from other stress conditions was the extensive induction of signaling and perception machinery playing a key role in effective resistance against pathogens. Interestingly, the signalling processes and immune response activated during BABA-IR in tomato differed from those in Arabidopsis with substantial enrichment of genes associated with jasmonic acid (JA) and ethylene (ET) signalling and no change in Asp levels. Our results revealed key differences between the effect of BABA on tomato and other model plants studied until now. Surprisingly, salicylic acid (SA) is not involved in BABA downstream signalization whereas ET and JA play a crucial role.

Cholesterol efflux pathways hinder KRAS-driven lung tumor progenitor cell expansion.

Cholesterol efflux pathways could be exploited in tumor biology to unravel cancer vulnerabilities. A mouse model of lung-tumor-bearing KRASG12D mutation with specific disruption of cholesterol efflux pathways in epithelial progenitor cells promoted tumor growth. Defective cholesterol efflux in epithelial progenitor cells governed their transcriptional landscape to support their expansion and create a pro-tolerogenic tumor microenvironment (TME). Overexpression of the apolipoprotein A-I, to raise HDL levels, protected these mice from tumor development and dire pathologic consequences. Mechanistically, HDL blunted a positive feedback loop between growth factor signaling pathways and cholesterol efflux pathways that cancer cells hijack to expand. Cholesterol removal therapy with cyclodextrin reduced tumor burden in progressing tumor by suppressing the proliferation and expansion of epithelial progenitor cells of tumor origin. Local and systemic perturbations of cholesterol efflux pathways were confirmed in human lung adenocarcinoma (LUAD). Our results position cholesterol removal therapy as a putative metabolic target in lung cancer progenitor cells.

The MIR34B/C genomic region contains multiple potential regulators of multiciliogenesis.

The MIR449 genomic locus encompasses several regulators of multiciliated cell (MCC) formation (multiciliogenesis). The miR-449 homologs miR-34b/c represent additional regulators of multiciliogenesis that are transcribed from another locus. Here, we characterized the expression of BTG4, LAYN, and HOATZ, located in the MIR34B/C locus using single-cell RNA-seq and super-resolution microscopy from human, mouse, or pig multiciliogenesis models. BTG4, LAYN, and HOATZ transcripts were expressed in both precursors and mature MCCs. The Layilin/LAYN protein was absent from primary cilia, but it was expressed in apical membrane regions or throughout motile cilia. LAYN silencing altered apical actin cap formation and multiciliogenesis. HOATZ protein was detected in primary cilia or throughout motile cilia. Altogether, our data suggest that the MIR34B/C locus may gather potential actors of multiciliogenesis.

The spatial landscape of gene expression isoforms in tissue sections.

In situ capturing technologies add tissue context to gene expression data, with the potential of providing a greater understanding of complex biological systems. However, splicing variants and full-length sequence heterogeneity cannot be characterized at spatial resolution with current transcriptome profiling methods. To that end, we introduce spatial isoform transcriptomics (SiT), an explorative method for characterizing spatial isoform variation and sequence heterogeneity using long-read sequencing. We show in mouse brain how SiT can be used to profile isoform expression and sequence heterogeneity in different areas of the tissue. SiT reveals regional isoform switching of Plp1 gene between different layers of the olfactory bulb, and the use of external single-cell data allows the nomination of cell types expressing each isoform. Furthermore, SiT identifies differential isoform usage for several major genes implicated in brain function (Snap25, Bin1, Gnas) that are independently validated by in situ sequencing. SiT also provides for the first time an in-depth A-to-I RNA editing map of the adult mouse brain. Data exploration can be performed through an online resource (https://www.isomics.eu), where isoform expression and RNA editing can be visualized in a spatial context.

Coupling live-cell imaging and in situ isolation of the same single cell to profile the transient states of predicted drug-tolerant cells.

Cell response variability is a starting point in cancer drug resistance that has been difficult to analyze because the tolerant cell states are short lived. Here, we present fate-seq, an approach to isolate single cells in their transient states of drug sensitivity or tolerance before profiling. The drug response is predicted in live cells, which are laser-captured by microdissection before any drug-induced change can alter their states. This framework enables the identification of the cell-state signatures causing differential cell decisions upon treatment. For complete details on the use and execution of this protocol, please refer to Meyer et al. (2020).

Unravelling the sex-specific diversity and functions of adrenal gland macrophages.

Despite the ubiquitous function of macrophages across the body, the diversity, origin, and function of adrenal gland macrophages remain largely unknown. We define the heterogeneity of adrenal gland immune cells using single-cell RNA sequencing and use genetic models to explore the developmental mechanisms yielding macrophage diversity. We define populations of monocyte-derived and embryonically seeded adrenal gland macrophages and identify a female-specific subset with low major histocompatibility complex (MHC) class II expression. In adulthood, monocyte recruitment dominates adrenal gland macrophage maintenance in female mice. Adrenal gland macrophage sub-tissular distribution follows a sex-dimorphic pattern, with MHC class IIlow macrophages located at the cortico-medullary junction. Macrophage sex dimorphism depends on the presence of the cortical X-zone. Adrenal gland macrophage depletion results in altered tissue homeostasis, modulated lipid metabolism, and decreased local aldosterone production during stress exposure. Overall, these data reveal the heterogeneity of adrenal gland macrophages and point toward sex-restricted distribution and functions of these cells.

The discovAIR project: a roadmap towards the Human Lung Cell Atlas.

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell-cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions.

TLR5 signalling is hyper-responsive in porcine cystic fibrosis airways epithelium.

Excessive lung inflammation and airway epithelium damage are hallmarks of cystic fibrosis (CF) disease. It is unclear whether lung inflammation is related to an intrinsic defect in the immune response or to chronic infection. We aimed to determine whether TLR5-mediated response is defective in the CF airway epithelium. We used a newborn CF pig model to study intrinsic alterations in CF airway epithelium innate immune response. Airway epithelial cells (AECs) were stimulated with flagellin or lipopolysaccharide to determine responses specific for TLR5 and TLR4, respectively. We observed a significant increase in cytokine secretion when CF AECs were stimulated with flagellin compared to wild type (WT) AECs. These results were recapitulated when AECs were treated with an inhibitor of CFTR channel activity. We show that TLR5-signalling is altered in CF lung epithelium at birth. Modulation of TLR5 signalling could contribute to better control the excessive inflammatory response observed in CF lungs.

Brown adipose tissue monocytes support tissue expansion.

Monocytes are part of the mononuclear phagocytic system. Monocytes play a central role during inflammatory conditions and a better understanding of their dynamics might open therapeutic opportunities. In the present study, we focused on the characterization and impact of monocytes on brown adipose tissue (BAT) functions during tissue remodeling. Single-cell RNA sequencing analysis of BAT immune cells uncovered a large diversity in monocyte and macrophage populations. Fate-mapping experiments demonstrated that the BAT macrophage pool requires constant replenishment from monocytes. Using a genetic model of BAT expansion, we found that brown fat monocyte numbers were selectively increased in this scenario. This observation was confirmed using a CCR2-binding radiotracer and positron emission tomography. Importantly, in line with their tissue recruitment, blood monocyte counts were decreased while bone marrow hematopoiesis was not affected. Monocyte depletion prevented brown adipose tissue expansion and altered its architecture. Podoplanin engagement is strictly required for BAT expansion. Together, these data redefine the diversity of immune cells in the BAT and emphasize the role of monocyte recruitment for tissue remodeling.

COVID-19 patients age, comorbidity profiles and clinical presentation related to the SARS-CoV-2 UK-variant spread in the Southeast of France.

The variant 20I/501Y.V1, associated to a higher risk of transmissibility, emerged in Nice city (Southeast of France, French Riviera) during January 2021. The pandemic has resumed late December 2020 in this area. A high incidence rate together with a fast turn-over of the main circulating variants, provided us the opportunity to analyze modifications in clinical profile and outcome traits. We performed an observational study in the University hospital of Nice from December 2020 to February 2021. We analyzed data of sequencing of SARS-CoV-2 from the sewage collector and PCR screening from all positive samples at the hospital. Then, we described the characteristics of all COVID-19 patients admitted in the emergency department (ED) (n = 1247) and those hospitalized in the infectious diseases ward or ICU (n = 232). The UK-variant was absent in this area in December, then increasingly spread in January representing 59% of the PCR screening performed mid-February. The rate of patients over 65 years admitted to the ED decreased from 63 to 50% (p = 0.001). The mean age of hospitalized patients in the infectious diseases ward decreased from 70.7 to 59.2 (p < 0.001) while the proportion of patients without comorbidity increased from 16 to 42% (p = 0.007). Spread of the UK-variant in the Southeast of France affects younger and healthier patients.

Monitoring SARS-CoV-2 variants alterations in Nice neighborhoods by wastewater nanopore sequencing.

BACKGROUND: Wastewater surveillance was proposed as an epidemiological tool to define the prevalence and evolution of the SARS-CoV-2 epidemics. However, most implemented SARS-CoV-2 wastewater surveillance projects were based on qPCR measurement of virus titers and did not address the mutational spectrum of SARS-CoV-2 circulating in the population. METHODS: We have implemented a nanopore RNA sequencing monitoring system in the city of Nice (France, 550,000 inhabitants). Between October 2020 and March 2021, we monthly analyzed the SARS-CoV-2 variants in 113 wastewater samples collected in the main wastewater treatment plant and 20 neighborhoods. FINDINGS: We initially detected the lineages predominant in Europe at the end of 2020 (B.1.160, B.1.177, B.1.367, B.1.474, and B.1.221). In January, a localized emergence of a variant (Spike:A522S) of the B.1.1.7 lineage occurred in one neighborhood. It rapidly spread and became dominant all over the city. Other variants of concern (B.1.351, P.1) were also detected in some neighborhoods, but at low frequency. Comparison with individual clinical samples collected during the same week showed that wastewater sequencing correctly identified the same lineages as those found in COVID-19 patients. INTERPRETATION: Wastewater sequencing allowed to document the diversity of SARS-CoV-2 sequences within the different neighborhoods of the city of Nice. Our results illustrate how sequencing of sewage samples can be used to track pathogen sequence diversity in the current pandemics and in future infectious disease outbreaks. TRANSLATION: For the French translation of the abstract see Supplementary Materials section.

Airway Administration of Flagellin Regulates the Inflammatory Response to Pseudomonas aeruginosa.

Excessive lung inflammation and airway epithelial damage are hallmarks of human inflammatory lung diseases, such as cystic fibrosis (CF). Enhancement of innate immunity provides protection against pathogens while reducing lung-damaging inflammation. However, the mechanisms underlying innate immunity-mediated protection in the lung remain mysterious, in part because of the lack of appropriate animal models for these human diseases. TLR5 (Toll-like receptor 5) stimulation by its specific ligand, the bacterial protein flagellin, has been proposed to enhance protection against several respiratory infectious diseases, although other cellular events, such as calcium signaling, may also control the intensity of the innate immune response. Here, we investigated the molecular events prompted by stimulation with flagellin and its role in regulating innate immunity in the lung of the pig, which is anatomically and genetically more similar to humans than rodent models. We found that flagellin treatment modulated NF-κB signaling and intracellular calcium homeostasis in airway epithelial cells. Flagellin pretreatment reduced the NF-κB nuclear translocation and the expression of proinflammatory cytokines to a second flagellin stimulus as well as to Pseudomonas aeruginosa infection. Moreover, in vivo administration of flagellin decreased the severity of P. aeruginosa-induced pneumonia. Then we confirmed these beneficial effects of flagellin in a pathological model of CF by using ex vivo precision-cut lung slices from a CF pigz model. These results provide evidence that flagellin treatment contributes to a better regulation of the inflammatory response in inflammatory lung diseases such as CF.

Cutaneous Squamous Cell Carcinoma Development Is Associated with a Temporal Infiltration of ILC1 and NK Cells with Immune Dysfunctions.

NK cells and tissue-resident innate lymphoid cells (ILCs) are innate effectors found in the skin. To investigate their temporal dynamics and specific functions throughout the development of cutaneous squamous cell carcinoma (cSCC), we combined transcriptomic and immunophenotyping analyses in mouse and human cSCCs. We identified an infiltration of NK cells and ILC1s as well as the presence of a few ILC3s. Adoptive transfer of NK cells in NK cell‒ and ILC-deficient Nfil3-/- mice revealed a role for NK cells in early control of cSCC. During tumor progression, we identified a population skewing with the infiltration of atypical ILC1 secreting inflammatory cytokines but reduced levels of IFN-γ at the papilloma stage. NK cells and ILC1s were functionally impaired, with reduced cytotoxicity and IFN-γ secretion associated with the downregulation of activating receptors. They also showed a high degree of heterogeneity in mouse and human cSCCs with the expression of several markers of exhaustion, including TIGIT on NK cells and PD-1 and TIM-3 on ILC1s. Our data show an enrichment in inflammatory ILC1 at the precancerous stage together with impaired antitumor functions in NK cells and ILC1 that could contribute to the development of cSCC and thus suggest that future immunotherapies should take both ILC populations into account.

Versatile and flexible microfluidic qPCR test for high-throughput SARS-CoV-2 and cellular response detection in nasopharyngeal swab samples.

The emergence and quick spread of SARS-CoV-2 has pointed at a low capacity response for testing large populations in many countries, in line of material, technical and staff limitations. The traditional RT-qPCR diagnostic test remains the reference method and is by far the most widely used test. These assays are limited to a few probe sets, require large sample PCR reaction volumes, along with an expensive and time-consuming RNA extraction step. Here we describe a quantitative nanofluidic assay that overcomes some of these shortcomings, based on the BiomarkTM instrument from Fluidigm. This system offers the possibility of performing 4608 qPCR end-points in a single run, equivalent to 192 clinical samples combined with 12 pairs of primers/probe sets in duplicate, thus allowing the monitoring of SARS-CoV-2 including the detection of specific SARS-CoV-2 variants, as well as the detection other pathogens and/or host cellular responses (virus receptors, response markers, microRNAs). The 10 nL-range volume of BiomarkTM reactions is compatible with sensitive and reproducible reactions that can be easily and cost-effectively adapted to various RT-qPCR configurations and sets of primers/probe. Finally, we also evaluated the use of inactivating lysis buffers composed of various detergents in the presence or absence of proteinase K to assess the compatibility of these buffers with a direct reverse transcription enzymatic step and we propose several protocols, bypassing the need for RNA purification. We advocate that the combined utilization of an optimized processing buffer and a high-throughput real-time PCR device would contribute to improve the turn-around-time to deliver the test results to patients and increase the SARS-CoV-2 testing capacities.

Correction: The nuclear hypoxia-regulated NLUCAT1 long non-coding RNA contributes to an aggressive phenotype in lung adenocarcinoma through regulation of oxidative stress.

Lung cancer is the leading cause of cancer death worldwide, with poor prognosis and a high rate of recurrence despite early surgical removal. Hypoxic regions within tumors represent sources of aggressiveness and resistance to therapy. Although long non-coding RNAs (lncRNAs) are increasingly recognized as major gene expression regulators, their regulation and function following hypoxic stress are still largely unexplored. Combining profiling studies on early-stage lung adenocarcinoma (LUAD) biopsies and on A549 LUAD cell lines cultured in normoxic or hypoxic conditions, we identified a subset of lncRNAs that are both correlated with the hypoxic status of tumors and regulated by hypoxia in vitro. We focused on a new transcript, Nuclear LUCAT1 (NLUCAT1), which is strongly upregulated by hypoxia in vitro and correlated with hypoxic markers and poor prognosis in LUADs. Full molecular characterization showed that NLUCAT1 is a large nuclear transcript composed of six exons and mainly regulated by NF-κB and NRF2 transcription factors. CRISPR-Cas9-mediated invalidation of NLUCAT1 revealed a decrease in proliferative and invasive properties, an increase in oxidative stress and a higher sensitivity to cisplatin-induced apoptosis. Transcriptome analysis of NLUCAT1-deficient cells showed repressed genes within the antioxidant and/or cisplatin-response networks. We demonstrated that the concomitant knockdown of four of these genes products, GPX2, GLRX, ALDH3A1, and PDK4, significantly increased ROS-dependent caspase activation, thus partially mimicking the consequences of NLUCAT1 inactivation in LUAD cells. Overall, we demonstrate that NLUCAT1 contributes to an aggressive phenotype in early-stage hypoxic tumors, suggesting it may represent a new potential therapeutic target in LUADs.

Single-cell RNA sequencing reveals intratumoral heterogeneity in primary uveal melanomas and identifies HES6 as a driver of the metastatic disease.

Intratumor heterogeneity has been recognized in numerous cancers as a major source of metastatic dissemination. In uveal melanomas, the existence and identity of specific subpopulations, their biological function and their contribution to metastasis remain unknown. Here, in multiscale analyses using single-cell RNA sequencing of six different primary uveal melanomas, we uncover an intratumoral heterogeneity at the genomic and transcriptomic level. We identify distinct transcriptional cell states and diverse tumor-associated populations in a subset of the samples. We also decipher a gene regulatory network underlying an invasive and poor prognosis state driven in part by the transcription factor HES6. HES6 heterogenous expression has been validated by RNAscope assays within primary human uveal melanomas, which further unveils the existence of these cells conveying a dismal prognosis in tumors diagnosed with a favorable outcome using bulk analyses. Depletion of HES6 impairs proliferation, migration and metastatic dissemination in vitro and in vivo using the chick chorioallantoic membrane assay, demonstrating the essential role of HES6 in uveal melanomas. Thus, single-cell analysis offers an unprecedented view of primary uveal melanoma heterogeneity, identifies bona fide biomarkers for metastatic cells in the primary tumor, and reveals targetable modules driving growth and metastasis formation. Significantly, our findings demonstrate that HES6 is a valid target to stop uveal melanoma progression.